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Abstract 

We present hydrokinetic Lattice Boltzmann and Molecular Dynamics simulations of capillary 
filling of high- wetting fluids in nano-channels, which provide clear evidence of the formation of 
thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor 
films is found to obey the Lucas- Washburn law as the main capillary front, z 2 (t) oc t, although 
with a larger prefactor, which we find to take the same value for both geometries under inspection. 
Both hydrokinetic and Molecular Dynamics approaches indicate a precursor film thickness of the 
order of one tenth of the capillary diameter. The quantitative agreement between the hydrokinetic 
and atomistic methods indicates that the formation and propagation of thin precursors can be 
handled at a mesoscopic/hydrokinetic level, thereby opening the possibility of using hydrokinetic 
methods to space-time scales and complex geometries of direct experimental relevance. 

PACS numbers: 
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Micro- and nanohydrodynamic flows are playing an emerging role for many applications 
in material science, chemistry and biology, and consequently they cast a pressing demand 
for a deeper and better understanding of the basic mechanisms of fluid flow at micro and 

nanoscales n a a a b, a a, as well as for the development of corresponding efficient 

computational tools. The formation of thin precursor films in capillary experiments with 
highly wetting fluids (near-zero contact angle) has been reported by a number of experiments 
and theoretical works ay , [h]] , mostly in connection with droplet spreading experiments. In 
spite of its conceptual and practical importance, however, much less work has been directed 
to the study of precursor films for the case of capillary filling. Indeed, this is of great interest 
because, by forming and propagating ahead of the main capillary front, thin-film precursors 
may manage to hide the chemical/geometrical details of the nanochannel walls, thereby 
exerting a major influence on the efficiency of nanochannel-coating strategies [ll|, [ijj]. In 
this realm, the continuum assumption behind the macroscopic description of fluid flow goes 
often under question, typical cases in point being slip-flow at solid walls and moving contact 
line of liquid/gas interface on solid walls [l|, [§]]. In order to keep a continuum description 
at nanoscopic scales and close to the boundaries, the hydrodynamic equations are usually 
enriched with generalised boundary conditions, designed in such a way as to collect the 
complex physics of fluid-wall interactions into a few effective parameters. A more radical 
approach is to quit the continuum level and turn directly to the atomistic description of fluid 
flows as a collection of moving molecules This approach is computationally demanding, 
thus preventing the attainment of space and time macroscopic scales of experimental interest. 
In recent years, a third, intermediate, alternative is rapidly materialising, in the form of 
minimal lattice versions of the Boltzmann kinetic equation [if], \v\ . By definition, such a 
mesoscopic approach is best suited to situations where molecular details, while sufficiently 
important to require substantial amendments of the continuum assumption, still possess a 
sufficient degree of universality to allow general continuum symmetries to survive, a situation 
that we shall dub supra-molecular for simplicity. Lacking a rigorous bottom-up derivation, 
the validity of the hydro-kinetic approach for supra-molecular physics must be assessed case- 
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by-case, a program which is already counting a number of recent successes 
The aim of this paper is to study another potentially 'supramolecular' situation, i.e. the 
formation and propagation of precursor films in capillary filling at nanoscopic scales. We 
provide quantitative evidences that the formation and propagation of a precursor film is 



3 



governed by kinetic-hydrodynamical laws, and that the emergent asymptotic dynamics for 
the head of the film follows a square root law -in time- as the bulk fluid. Nevertheless, the 
prefactors of the two behaviours are different, eventually leading to a stable infinitely long 
thin precursor film ahead of the front. We employ both MD and hydrokinetic simulations, 
finding quantitative agreement for both bulk quantities and local density profiles, at all 
times during the capillary filling process. This lends further support to the idea that many 
supramolecular problems are efficiently handled with mesoscopic/hydrokinetic methods. 



In this work, we use the multicomponent LB model proposed by Shan and Chen [23 ]. 
This model caters for an arbitrary number of components, with different molecular mass 
and reads as follows: 

ft (x + *M,t + At) - /*(*,*) = -—[//»& t) - f* {eq \x,t)] (1) 

where f^(x,t) is the probability density function associated with a mesoscopic particle of 
velocity for the kth component, r k is a mean collision time of the kth component (At being 
the time-step associated with particle free-streaming), and (x,t) is the corresponding 
Maxwell equilibrium function. The interaction force between particles is the sum of a bulk 
and a wall components: 

s 

F bk (x) = -n k (x)^^G kk n k (x)(x + Ax); F wk (x) = -n k (x) 2J g kw n w {x){x + Ax) (2) 

x fe=l X 

where n k is the local density and the interaction matrix G kk is symmetric . In our model, 
G kk = g kk , (for \x' — x\ = 1);, G kk = g kk /^-, (for \x' — x\ = a/2; and otherwise. At the 
fluid/solid interface, the wall is regarded as a (solid) phase with a constant number density 
n w and g kw is the interaction strength between component k and the wall. By adjusting the 



parameters g kw and n w , different wettability can be obtained [251 ]. This approach allows 
the definition of a static contact angle 8, spanning the entire range 9 e [0° : 180°] [26]. To 
the purpose of analysing the physics of film precursors, it is important to notice that , in 
the limit where the hard-core repulsion is negligible, both the Shan-Chen pseudo-potential 
and van der Waals interactions predict a non-ideal equation of state, in which the leading 



correction to the ideal pressure is oc p 2 . Concerning the MD simulation 22|, see fig. [H, 
fluid atoms interact through a LJ potential, Ulj{t) = 4e[(a/r) 12 — (cr/r) 6 ], where e = 1.4 
and a = 1.0. The capillary walls are represented by atoms forming a triangular lattice with 
spacing 1.0 in units of the liquid atom diameter a. The wall atoms may fluctuate around 



their equilibrium positions at R+a, subject to a finitely extensible non-linear elastic (FENE) 
potential Ufene = l5e w Rlln(l — r 2 /Rl), i? = 1-5. Here e w = l.OksT, where ks denotes 
the Boltzmann constant and T is the temperature of the system. In addition, the wall 
atoms interact through a LJ potential, Ulj{t) = 4:6 ww [(a ww /r) 12 ((a ww /r) 6 ], where e ww = 1.0 
and a ww = 0.8. This choice guarantees no penetration of liquid particles through the wall. 
Molecules are advanced in time via the velocity- Verlet algorithm and a dissipative-particle- 



dynamics thermostat, with friction parameter £ = 0.5, and a thermostat cutoff r c = 2.5cr [15]. 
The integration time step is 8t = O.Olto, where to = ^/ma 2 / A8ksT = l/y/48 is the basic 
time unit and we have taken the particle mass m = 1 and ksT = 1. 

We consider a capillary filling experiment, whereby a dense fluid, of density and dynamic 
viscosity pi,/ii, penetrates into a channel filled up by a lighter fluid, p2,A i 2, see fig. [TJ For 



this kind o 
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fluid flow, the Lucas- Washburn law is expected to hold, at least at macroscopic 
I4J] and in the limit [i\ ^> /12. Recently, the same law has been observed even 
in nanoscopic experiments 24|. In these limits, the law governing the position, z(t), of the 
macroscopic meniscus is: z (t) - z 2 (0) = lE0fil tj W here 7 is the surface tension between 
liquid and gas, 9 is the static contact angle, \x is the liquid viscosity, H is the channel height 
and the factor C depends on the flow geometry (in the present geometry Clb = 3 ; Cmd = 
2). The geometry we are going to investigate is depicted in fig. for both models. It is 
important to underline that in the LB case, we simulate two immiscible fluids, without any 
phase transition. 

Since binary LB methods do not easily support high density ratios between the two 
species, we impose the correct ratio between the two dynamics viscosities through an appro- 
priate choice of the kinematic viscosities. The chosen parameters correspond to an averaged 
capillary number Ca w 3 10~ 2 and Ca ~ 0.1 for LB and MD respectively. 

In order to expose the universal character of the phenomenon, results are presented in 
natural units, namely, space is measured in units of the capillary size, l cap and time in units of 
the capillary time t cap = l C a P /Vca P , where V cap = 7/// is the capillary speed and l cap = H/Clb 
for LB and l cap = R/Cmd for MD. In these units, the Lucas- Washburn law takes a very 
simple universal form 
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z 2 +t, (3) 



where we have inserted the value of cos(8) = 1, corresponding to complete wetting. As to 
the bulk front position, figl2^ shows that both MD and LB results superpose with the law 
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(a) Reservoir 

FIG. 1: Sketch of the geometry used for the description of the capillary imbibition in the LB and MD 
simulations, (a) The 2 dimensional geometry, with length L\ + L2 and width H, is divided in two 
parts. The left part has top and bottom periodic boundary conditions, so as to support a perfectly 
flat gas- liquid interface, mimicking an "infinite reservoir". In the right half, of length L%, there is 
the actual capillary: the top and bottom boundary conditions are solid wall, with a given contact 
angle 6. Periodic boundary conditions are also imposed at the inlet and outlet sides. The main LB 
parameters are: H = ny = 40, Li = nz = 170, pi = l;p2 = 0.35, /^i = 0.66,^2 = 0.014,7 = 0.016 
where H is the channel height, L2 is the channel length, p2 and p\ the gas and liquid densities 
respectively: p^, k = 1, 2 the dynamic viscosities, and 7 the surface tension, (b) Snapshot of fluid 
imbibition for MD in the capillary at time t = 1300 MD time-steps. The fluid is in equilibrium with 
its vapour. Fluid atoms are in blue. Vapour is yellow, tube walls are red and the precursor is green. 
One distinguishes between vapour and precursor, subject to the radial distance of the respective 
atoms from the tube wall, if a certa in p article has no contact with the wall, it is deemed 'vapour'. 



The MD parameters are as follows 



22J: R = 11a, L = 55a, p x = 0.774, p = 6.3, 7 = 0.735, a = 1, 



where R is the capillary radius and L its lenght. 



d2J), while the precursor position develops a faster dynamics, fitted by the relation: 
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Similar speed-up of the precursor has been reported also in different experimental situations 
[lo| . The precursor is here denned through the density profile, ~p(z), averaged over the 
direction across the channel. In figure EJo, we show p(z) at various instants (always in 
capillary units), for both MD and LB simulations. The relatively high average density ~p(z) 
near the wall, at later times, witnesses the presence of a precursor film attached to the wall. 
From this figure, it is appreciated that quantitative agreement between MD and LB is found 
also at the level of the spatial profile of the density field. This is plausible, since the LB 
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FIG. 2: Dynamics of the bulk and precursor meniscus, (a) Position of the liquid meniscus z (t) 
for LB and MD simulations. The position of the precursor film, Zp rec (t) is also plotted for both 
models. z prec is defined as the rightmost location with density p = Pbuik/3- All quantities are given 
in natural "capillary" units (see text). The asymptotic (t > 15i cap ) rise of both precursor and 
bulk menisci follows a Lucas- Washburn law, with different prefactors (see the two straight lines). 
Notably, the precursor film is found to proceed with the l&w.Zp rec (t) = 1.35t. (b) Profiles of the 
average fluid density ~p(z) in the capillary at various times for LB and MD models. 

simulations operate on similar principles as the MD ones, namely the fluid-wall interactions 
are not expressed in terms of boundary conditions on the contact angle, like in continuum 
methods, but rather in terms of fluid-solid (pseudo)-potentials. In particular, the degree of 
hydrophob/philicity of the two species can be tuned independently, by choosing different 
values of the fluid-solid interaction strengths (for details see 221]). In figure Eb, we show the 
position of the advancing front and the precursor film as a function of time, for both MD 
and LB simulations. Even though the averaged capillary numbers are not exactly the same, 
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(a) 

FIG. 3: Structure of the precursor film, (a) Snapshot of the interface profile for LB and MD models 
at time t = 20. Both LB and MD profiles are obtained as isolines at Pbuik/p = 3. We also show 
the circular fit (bold line) , valid for the menisci in the bulk region and the 1/z profile as suggested 
by lubrication theory. Notice the clear departure of the interface form a circular fit, due to the 
different physics close to the wall. Natural "capillary" units are used (see text), (b): The interface 
profile at a later stage of capillary filling (symbol: square). The circle-straight line corresponds to 
the circular fit of the bulk (0.15 < y < 1.5). Notice the long extension of the precursor, here defined 
as the spatial interval where the profile departs from the circular fit. The length of the precursor 
film is consistent with the relation L p ~ Ca~ 1 , as found in a recent nanoscale experiment [9(] (not 
shown). In the inset, we plot the value of the macroscopic angle defined by the interface curvature 
at varying the distance from the wall. 



a pretty good agreement between LB and MD is observed. In particular, in both cases, 
the precursor is also found to obey a \ft scaling in time, although with a larger prefactor 
than the front. As a result, the relative distance between the two keeps growing in time, 
with the precursors serving as a sort of 'carpet', hiding the chemical structure of the wall 



to the advancing front. As recently pointed out [22|, |24(, MD simulations and experiments 
indicate that the Lucas- Washburn law persists down to nanoscopic scales. The fact that 
the precursor dynamics quantitatively matches with mesoscopic simulations, suggests that 
the precursor physics also shows the same kind of nanoscopic persistence. Concerning the 
instantaneous interface profile, in Fig. [3^ we show the structure of the precursor by analysing 
fluid menisci cuts for different density levels. The agreement between the two methods is 
again quantitative. In LB, the interface turns out to be about 5 Az (leading to the rough 
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estimate, Az ~ 0.2nm, considering an interface thickness for water/air at room temperature 
A£ ^ Inm). The film profile, defined by the isoline at density p = pbuik/3, splits in two 
regions. In the first one, reaching up to y ~ 7Az, the front is fitted by the expected bulk 
circle profile. In the second one, extending from y w 7Az to the wall, the film profile is 
fitted by the analytical expression ~ Ca _1 a 2 /z, where a is a characteristic molecular size, 
which in our case is set equal to a = 1. The analytical profile was proposed for the case 
of "maximal film" (perfect wetting) under the lubrication approximation, with Van der 
Waals interactions between fluid and walls. Quantitative agreement is again found between 
LB, MD and the theoretical treatment. 

In order to further analyse the previous results, we study the asymptotic case, i.e. a 
situation with a very small fluid velocity (quasi-static), obtained in a longer channel (L = 
nz = 1500A2;) with the same geometry depicted in fig. [TJ To be noted that such domain 
is too large for an atomistic method and, therefore, we have simulated it with LB only. In 
such a case, the departure between the bulk profile and the precursor shape is even more 
pronounced, due to the longer elongation of the film close to the wall. In the inset of the 
same figure, we also show the macroscopic angle formed by the interface as a function of 
the distance from the wall. The bulk interface is fitted by a circle extrapolating to the 
wall with a contact angle 8 = 0. It is remarkable that the presence of the precursor film 
allows the apparent contact angle to vanish less sharply than in the ideal circular case, as 
visible in the inset of fig. [3b. The deviation from the circular fit can not be due to a simple 
dynamical distortion of the apparent contact angle, since this would lead to an increase of 
the macroscopic angle. Therefore, this deviation can be interpreted as a signature of the 
onset of precursor physics. 

Summarising, we have reported quantitative evidence of the formation and the dynamics 
of precursor films in capillary filling with highly wettable boundaries. The precursor shape 
shows persistent deviation from an ideal circular meniscus, due to the nanoscopic distortion 
induced by the interactions with the walls. This is likely connected to the disjoining pressure 
induced by van der Waals interactions between fluid and solid [2]. The precursor follows 
the same Lucas- Washburn functional law as the bulk meniscus, although with a different 
prefactor, see eq. (Jl|). This prefactor has been found to be the same for both geometries 
investigated, while its universality with respect to other physical parameters remains to be 
explored. Our findings are supported by direct comparison between LB, MD simulations 
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and theoretical predictions. On the practical side, the hydrokinetic method presented in 
this work offers the opportunity to perform very efficient numerical simulations of precursor 
formation and propagation in complex geometries of direct interest for the optimal design of 
micro- and nanofluidic devices. Furthermore, a very intringuing result is that the precursor 
film near-to-the-wall is found to cause a distorsion of the apparent contact angle, which 
suggests the possibility to design devices with increased efficiency in presence of roughness. 
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